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Abstract 

The prospect of base station (BS) cooperation leading to joint combining at widely separated antennas has 
led to increased interest in macrodiversity systems, where both sources and receive antennas are geographically 
distributed. In this scenario, little is known analytically about channel capacity since the channel matrices have a 
very general form where each path may have a different power. Hence, in this paper we consider the ergodic sum 
capacity of a macrodiversity MIMO system with arbitrary numbers of sources and receive antennas operating over 
Rayleigh fading channels. For this system, we compute the exact ergodic capacity for a two-source system and a 
compact approximation for the general system, which is shown to be very accurate over a wide range of cases. 
Finally, we develop a highly simplified upper-bound which leads to insights into the relationship between capacity 
and the channel powers. Results are verified by Monte Carlo simulations and the impact on capacity of various 
channel power profiles is investigated. 

Index Terms 

Macrodiversity, MIMO, MIMO-MAC, Capacity, Sum-rate, Network MIMO, CoMR DAS, Rayleigh fading. 

I. INTRODUCTION 

With the advent of network multiple input multiple output (MIMO) [1], base station (BS) collaboration 
[0 and cooperative MIMO (3), it is becoming more common to consider MIMO links where the receive 
array, transmit array or both are widely separated. In these scenarios, individual antennas from a single 
effective array may be separated by a considerable distance. When both transmitter and receiver have 
distributed antennas, we refer to the link as a macrodiversity MIMO link. Little is known analytically 
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about such links, despite their growing importance in research - flVj and standards where coordinated 
multipoint transmission (CoMP) is part of 3GPP LTE Advanced. 

Some analytical progress in this area has been made recently in the performance analysis of linear 
combining for macrodiversity systems in Rayleigh fading [[8]|, 0. However, there appears to be no 
work currently available on the capacity of general systems of this type. Similar work includes the 
capacity analysis of Rayleigh channels with a two-sided Kronecker correlation structure ifTOl . However, the 
Kronecker structure is much too restrictive for a macrodiversity layout and such results cannot be leveraged 
here. Also, there is interesting work on system capacity for particular cellular structures, including Wyner's 
circular cellular array model [6] and the infinite linear cell-array model [7]. Despite these contributions, 
the general macrodiversity model appears difficult to handle. The analytical difficulties are caused by the 
geographical separation of the antennas which results in different entries of the channel matrix having 
different powers with an arbitrary pattern. Also, these powers can vary enormously when shadowing and 
path loss are considered. Note that this type of channel model also occurs in the work of IfTTTl . 
In this paper, we consider a macrodiversity MIMO multiple access channel (MIMO-MAC) where all 
sources and receive antennas are widely separated and all links experience independent Rayleigh fading. 
For this system, we consider the ergodic sum capacity, under the assumption of no channel state infor- 
mation (CSI) at the transmitters. For two sources, we derive the exact ergodic sum capacity. The result is 
given in closed form, but the details are complicated and for more than two sources, it would appear that 
an exact approach is too complex to be useful. Hence, we develop an approximation and a bound for the 
general case. The first technique is very accurate, but the functional form is awkward to interpret. Hence, 
a second, less accurate but simple bound is developed which has a familiar and appealing structure. This 
bound leads to insight into capacity behavior and its relationship with the channel powers. In |[T2ll . we 
presented a preliminary study of this problem, which focussed on the approximation for the general case. 
In this paper, we have extended the conference version to include the exact two source results, correlated 
channels, full mathematical details (see Sec. Hill) , a motivation for the approximate analysis (see Appendix 
ID1) and a much wider range of scenarios, power profiles and discussion in the results section. 
Note that, the methodology developed is for the case of arbitrary powers for the entries in the channel 
matrix. There is no restriction due to particular cellular structures. Hence, the results and techniques may 
also have applications in multivariate statistics. 

The rest of the paper is laid out as follows. Section [TT] describes the system model and Sec. [TIT] gives some 



mathematical preliminaries required in the analysis. Section [IV] provides an exact analysis for the case of 
two source antennas. Sections [V] and [VI] consider the case of arbitrary numbers of sources and develop 
accurate approximations and bounds on capacity. Results and conclusions appear in Sees. IVIII and [Villi 

II. System Model 

Consider a MIMO-MAC link with M base stations and W users operating over a Rayleigh channel 
where BS i has n R . receive antennas and user i has rij antennas. The total number of receive antennas is 
denoted tir = YliLi n Ri an( ^ me tota l number of transmit antennas is denoted N = $^ i=1 rii. An example 
of such a system is shown in Fig. [TJ where three BSs are linked by a backhaul processing unit (BPU) 
and communicate with multiple, mobile users. All channels are considered to be independent since the 
correlated channel scenario can be transformed into the independent case as shown in Sec. III-AI The 
system equation is given by 

r = Hs + n, (1) 

where r is the C nRXl receive vector, s is the combined C Nxl transmitted vector from the W users, n is an 
additive white Gaussian noise vector, n ~ CM (0, a 2 1), and H E C UrxN is the composite channel matrix 
containing the W channel matrices from the W users. The ergodic sum capacity of the link depends on 
the availability of channel state information (CSI) at the transmitter side. In particular, if no CSI at the 
transmitter is assumed, the corresponding ergodic sum capacity is fl3] pp. 57] 



E{C} = E\\og 2 



I + - n HH H 



a 2 



(2) 



where i?{|sj| 2 } = 1, i = 1,2, . . . , N, is the power of each transmitted symbol. It is convenient to label 
each column of H as hi, i = 1, 2, . . . , N, so that H = (hi,h 2 , . . . ,hjy). The covariance matrix of h k 

is defined by P k = E ^hkhfj and Pk = diag (Pik, P2k, ■ ■ ■ , Pn R k)- Hence, the ik th element of H is 

i 

CM (0, Pik). Using this notation, we can also express hk as hk = Pk u k, where Uk ~ CM (0, J). Note that, 
for convenience, all the power information is contained in the Pk matrices so that there is no normalization 
of the channel and, in ©, the scaling factor in the capacity equation is simply I /a 2 . 



A. Correlated Channels 

Consider the general scenario where sources and/or BSs have multiple co-located antennas for trans- 
mission and reception. Here, spatial correlation may be present due to the co-located antennas |[T3l . [[141 . 




Fig. 1. A network MIMO system with a 3 sector cluster. To reduce the clutter, only paths from a single source are shown. 



H 



(3) 



If a Kronecker correlation model is assumed, then the composite channel matrix is given by 

f Ri \/H w>u ... H wiw \/Ri \ 

R? M J - H WiMW J \ o Rf w J 

where the C nfl fc xni matrix, H W:ik , has iid elements since all the channel powers from user k to BS i are 
the same. The matrix R ri is the receive correlation matrix at BS i and the matrix R t k is the transmit 
correlation matrix at source k as defined in |fl3l . Using the spectral decompositions, R„ = $ r iA ri $^ and 
Rtk = ^tk^-tk^u: an d substituting © into © it is easily shown that the capacity with the channel in © 
is statistically identical to the capacity with channel 

'Ai \{H w ,u ... H w , iw \fkl 

w .All 

H 



H 



(4) 







A| 



Denoting © by = H w Af , we see that correlation is equivalent to a scaling of the channel by the 
relevant eigenvalues in A r and A t . In particular, the (u, v) th element of H has power A rtUU A t , vv P uv , where 
P uv is the single link power from transmit antenna v to receive antenna u. Hence, correlation can be 
handled by the same methodology developed in Sees. ITVWll with suitably scaled power values Q. 



III. Preliminaries 

In this section we derive some useful results which will be used extensively throughout the paper. 

'Arbitrary fixed transmit power control techniques can also be handled in the same way as for the correlated scenario. 
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Lemma 1. Let X be an n x n complex random matrix with, 

( 



A=E{XoX) 



£{|X n | 2 } E{\X 12 \ 2 } . 
E{\X 2l \ 2 } E{\X 22 \ 2 } . 



(5) 



\E{\X nl \ 2 } E{\X n2 \ 2 } ...J 
where o represents the Hadamard product. With this notation, the following identity holds. 



E{\X H X\} = perm (A), 

where perm(.) is the permanent of a square matrix defined in 4731/ . 

Proof: From the definition of the determinant of a generic matrix, X = {^i,fc} ifc=1 n , we have 



(6) 



E{\X H X\] = E 



5^sgn(^)II X /^ 



i=l 



(7) 



where a = (eri, a 2 , ■ ■ ■ , cr n ) is a permutation of the integers 1, . . . , n, the sum is over all permutations, 
and sgn(a) denotes the sign of the permutation. The permutation, /i, in the second summation of £7J 
is defined similarly. Since all the elements of X are independent, the only terms giving non-zero value 
expectations are niLi X au i Yii=i ^h,h where permutation a = b. Hence, using the permanent definition 
in |fT5l we have 



E{\X H X\} = J2U A ^=V Qrm W- 



(8) 



a i=l 



Corollary 1. Let X be an m x n random matrix with, E {X o X} = A, where A is an mxn deterministic 
matrix and m > n. Then, the following identity holds. 



E{ \X"X\ \ = Perm (A) 



where Perm (A) is the permanent of the rectangular matrix A as defined in 4751/ . 



(9) 



Proof: Using the Cauchy-Binet formula for the determinant of the product of two rectangular matrices, 
we can expand IX^Xl as a sum of products of two square matrices. Each product of square matrices 
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can be evaluated using Lemma [TJ The resulting expression is seen to be the permanent of the rectangular 
matrix, A, which completes the proof. ■ 

Corollary 2. Let X be an m x n random matrix with, E {X o X} = A, where A is an mxn deterministic 
matrix and m > n. If the m x m deterministic matrix E is diagonal, then the following identity holds. 

E { \X H EX\ } = Perm (HA) . (10) 



Proof: The result follows directly from Lemma [Q and the fact that S 2 o E 2 = E for any diagonal 
matrix. ■ 
Next, we give a definition for the elementary symmetric function (esf) of degree k in n variables, 
X 1 , X 2 , . . . ,X n |fT6l . Let e k X 2 , . . . , X n ) be the k th degree esf, then 



e/c (Xi , X 2 , . . . , X, t 



l<ii<Z 2 <— <ifc<n 



x h . . . x h . 



(11) 



It is apparent from (fTTT) that e X 2 , . . . , X n ) = 1 and e n X 2 , . . . , X n ) = X\X 2 . . . X n . In general, 
the esf of degree k in n variables, for any k < n, is formed by adding together all distinct products of k 
distinct variables. 

Lemma 2. / Tiol/ let X be an nxn complex symmetric positive definite matrix with eigenvalues Ai, . . . , A n . 
r/?en, following identity holds. 



e k (Xi, X 2 , . . . , \ n ) = Tr k (X) 



(12) 



where 



Tr fc (X) 



1 





1 < k < n 
k = 
k > n, 



(13) 



where a k n is an ordered subset of {n} — {1, . . . , n} of length k and the summation is over all such subsets. 
X rr k n denotes the principal submatrix of X formed by taking only the rows and columns indexed by 

Cfc.n- 



In general, X^' n denotes the submatrix of X formed by taking only the rows and columns indexed 
by a it n and fi£^ n respectively, where cr^ n and \i^ n are length £ subsets of {1,2, ...,n}. If either o^ n 
or ii^ n contains the complete set (i.e., a^ n = {1,2, ... ,n) or ^ n = {1,2, .. . , n}), the corresponding 
subscript/superscript may be dropped. When o> n = fi e n , only one subscript/superscript may be shown for 
brevity. 

Next, we present three axiomatic identities for permanents which are required in Sec. |V] 
• Axiom 1: Let A be an arbitrary m x n matrix, then 



^Perm^)^) = ^Perm ((A) aom ) = 1. 

a cr 

Axiom 2: Let A be an arbitrary m x n matrix, then 

£ Perm ((A)^) = £ Perm ((A'f k -" ) . 

a cr 

Axiom 3: For an empty matrix, A. 

Perm (A) = 1. 



(14) 



(15) 



(16) 



IV. Exact Small System Analysis 

In this section, we derive the exact ergodic sum capacity in © for the N = 2 case. This corresponds to 
two single antenna users or a single user with two distributed antennas. Here, the channel matrix becomes 
H = (hi,h 2 ) and it is straightforward to write © as 



E{C}\n2 = £ jln J + -^Vif J 



+ E I In 



(7 Z 



7+1 ( J+i-Mf ) Ma 



-i 

H \ i iff 



Ci + c 2 . 



(17) 



Both Ci and C 2 can be expressed as scalars ifTTll . |[T8l pp. 48], so the capacity analysis simply requires 

1 



d = E I In ( 1 + — /if/l! 



0" 1 



(18) 
(19) 



In order to facilitate our analysis, it is useful to avoid the logarithm representations in (TT8T ) and (fT9l ). We 
exchange logarithms for exponentials as follows. First, we note the identity, 



1 

a 



e~ at dt, for a > 0. (20) 



o 



Now equation (1201) can be used to find In a as below: 

d In a 



oo 

at 



e~ at dt, (21) 



o 



da 

In a pa />oo 

d\na= / e~ at dtda, (22) 
o Ji Jo 

f°° e~ l — e~ at 

lna= / dt. (23) 

Jo t 

This manipulation is useful because there are many results which can be applied to exponentials of 
quadratic forms, whereas few results exist for logarithms. As an example, using (1231) in (TT8T) gives 




00 e~ l - e 



-dt) . (24) 



t 

Note that a = 1 + -^h^hi has been used in (|23l) . Sice a > 1, it follows that the integrand in (|24|) is 
non-negative. Also, the expected value, C\, is clearly finite and so, by Fubini's theorem, the order of 
expectation and integration in d24l) may be interchanged. Using the Gaussian integral identity 0, the 
expectation in (|24|) can be computed to give 

Ci =i -r-im dt ' <25) 

where Ei = J + \P\. Hence, the log-exponential conversion in (|23l leads to a manageable integral for 
C\. Using the same approach and applying d23l in (fl9l) gives 

-*-^f (l+jzhihf} h 2 
C'i E{ I —-- ^ dt>- (26) 




The expectation in (|26l) has to be calculated in two stages. First, the expectation over h 2 can be solved 
using the Gaussian integral identity [9] and, with some simplifications, we arrive at 



oo e -t 



e -t f a 2 + h H hi 

C 2 = I —~E h { ± '-^ } dt, (27) 



where S 2 = I+-^P 2 - Interchange of the expectation and integral in (1271 ) follows from the same arguments 
used for G\. Equation (ITTT ) can be further simplified to give 



C 2 = ^cft - £ ft < -r / , — L -^dt \ . (28) 



ct2 -/o £|E 2 | (a 2 + /if E 2 % 
Defining the third term in (|28l) as it, the ergodic sum capacity, i£ (C) = C\ + C 2 , becomes 



{C} = ^{t U-b}, (29) 



oo r e -t e -t 



E 

where 

'»>=X [T-m ]dt <30) 

Substituting for S fc in d30l) and expanding (t |Sfc|) 1 gives 
where 



^ifc = Ffiii / p p-^T- (32) 



Note that the first e^/t term in (1301) cancels out with one of the terms in the partial fraction expansion 
leaving only the linear terms shown in the denominator of (TJTI) . The integrals in (|3TT) can be solved in 
closed form |fT9l to give 



1=1 



J^^^e^ f— J. (33) 



In order to compute If, we use flU Lemma 1] to give 



OO /"OO 



</0 



|S 2 | 00i 



d0 2 dfc, (34) 

01=0 



10 



where and z 2 = a 2 + hf S 2 1 hi. The expectation in (|34l ) can be solved as in d9]|, and 

with some manipulations we arrive at 



OO /*00 Q 




|j + tp 2 + e 1 p l p 2 + e 2 Pi 



d6 2 dt. 



(35) 



01=0 



In Appendix [7a1 If, in (1351 ) is calculated in closed form and the final result is given by 



1 i^^^iiki(M hm -N bik 
|PiP 2 l l 



J,: 



(36) 



I i=l ky^i l^i,k 

where M& ifci , N b . kl , J { and £ iM are given in (|94j (1951) , (1721) and (|82l respectively. Then, the final result 



becomes 



E{C} 



■y ( 2 nR 2 / 2 



|PlP 2 



"■fl "j? n R 
i=l fc^i Z^i,fc 



dki [M b . kl - N b . k 



Ji 



(37) 



V. Approximate General Analysis 

In this section, we present an approximate ergodic sum rate capacity analysis for the case where 

n-R > N > 2. Extending this to iV > n R is a simple extension of the current analysis. We use the 
following notation for the channel matrix, 

H=(H N ,h N ) (38a) 

= (H N „ u h N _ u h N J (38b) 

Hi. hi, h\ ..h.x) (38c) 



(hi,h 2 , . . . , Hn) 



(38d) 



where the n R x (k — 1) matrix, H k , comprises the k — 1 columns to the left of hk in H. Using the same 
representation as in (fTTl) . the ergodic sum capacity is defined by [|3l as, 



N 



E{C}ln2±J2 C >" 



(39) 



fc=i 
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where 



C k = El\n 



I rX(l+-^H k H" ) h k h» 



a 2 \ a 2 



-i 



H 



(40) 
(41) 



Applying ([23]) to gB gives 



C k 



oo -t 

E 

t 





a 2 I + H k H k 









dt. 



(42) 



where = I + In order to calculate the second term in (l42l) . the following expectation needs to 

be calculated, 



/*(*) 



1 



IE* 



a 2 I + H k H k 



a 2 I + H'X ] H k 



(43) 



Exact analysis of J fc (t) is cumbersome, and even the N = 2 case (see the calculation in (|35l) ) is 
complicated. Hence, we employ a Laplace type approximation |20], so that I k (t) can be approximated 
by 



h (t) 



1 E 


{ 


a 2 I + H k H k 


} 






r 




} 



(44) 



Note that the Laplace approximation is better known for ratios of scalar quadratic forms 11201 . However, 
the identity in both the numerator and denominator of (1431) can be expressed as the limit of a Wishart 
matrix as in [|22l . This gives (|43T) as ratio of determinants of matrix quadratic forms which in turn can be 
decomposed to give a product of scalar quadratic forms as in Appendix iDl and [23 J. Hence, the Laplace 
approximation for (1431) has some motivation in the work of |{20l . It can also be thought of as a first order 
delta expansion ll24l . From Appendix [Bj the expectation in the numerator of (l44l) is given by 



k— 1 

E { | a*I + H H k H k | } = E E Perm ((^r^ 1 ) M ^ 



(45) 



i=0 a 



where Q k is defined in (I100I ). From Appendix the expectation in the denominator of (1441) is given by 



|Efe| E 



riR 
1=0 



(46) 
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where (p ki is given in (II 101) . Therefore, I k (t) becomes 



hin*JH&- (47) 
e (Q k ) 



e(Q fc ) 

nr=i (* + ' 



(48) 



(49) 



where 

fe-i 



1=0 cr 

Note that cjjy > for all l,k from Descartes' rule of signs as all the coefficients of the monomial in the 
denominator of (l48l are positive. Also note that, from (II 101) . we have 6 (Q k ) = ip k Q. Applying (|49l ) in 
(1421) gives 

Using a partial fraction expansion for the product in the denominator of the second term of (15TI) gives 



"R 



t nr=i (* + u ki) t t + u k i 



where 



,feo 



1 _ fkna 

TfTlR 



(53) 



and 



ski 



Applying (|52j) in (02J) gives 



n R poo > 

Ck ^J£MLJ2 -^-dt (55) 



(54) 



</?fco 



i=i 

IIR 



J2Ckie 0Jkl E 1 (u kl ). (56) 
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Then, applying (|56l ) in (1391) gives the final approximate ergodic sum rate capacity as 



, N / n R \ 

*<^EfeE^* <<■*>)■ 



Note the simplicity of the general approximation in (1571) in comparison to the two-user exact results in 
073. 



VI. A Simple Capacity Bound 

In this section, we derive an extremely simple upper-bound for the ergodic capacity in ©. This provides 
a simpler relationship between the average link powers and ergodic sum capacity at the expense of a loss 
in accuracy. Using Jensen's inequality [1211 and 7 = \, equation © leads to 

E {C} < log 2 (E{\I + ^H H H\ }) . (58) 

From Appendix El (|58l) can be given as 

e {c} < i og2 (j2 E Perm ( p<T,JV ) ? ) > (59) 

\i=0 a ) 

= log 2 fe^J • (60) 

where P = (P^). The simplicity of (l59~l) is hidden by the permanent form. For small systems, expanding 
the permanent reveals the simple relationship between the upper bound and the channel powers. For 
n R = N = 2 and n R = N = 3, (|60l ) gives the upper bounds in (loTT ) and (|62l ) respectively. These bounds 
show the simple pattern where cross products of L channel powers scale the 7 L term. Hence, at low 
SNR where the 7 term is dominant, log 2 (1 + Prr), where Pp = ^ P^, is an approximation to (1601 ). 
Similarly, at high SNR, the 7^ term is dominant and log 2 (l + Perm (P) 7^) is an approximation. These 
approximations show that capacity is affected by the sum of the channel powers at low SNR, whereas at 
high SNR, the cross products of N powers becomes important. 

VII. Numerical and Simulation Results 

For the numerical results, we consider three distributed BSs with either a single receive antenna or 
two antennas. For a two-source system, we parameterize the system by three parameters (HI, AH. The 
average received signal to noise ratio is defined by p = Pr/cr 2 . In particular for a two-source system, 
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E {C} < log 2 (1 + 7 (Pn + P 12 + P 21 + P 22 ) + f (P11P22 + P12P21)) ■ (61) 

E {C} < log 2 (I + 7 (P n + P 12 + P 13 + P 21 + P 22 + P 23 + P31 + P 32 + P33) 

+ 7 2 (P11P22 + P11P52 + P21P12 + P21P32 + P31P12 + P31P22 + P11P23 + P11P33 + P21P13 

(62) 

+ P21P33 + P31P13 + P31P23 + P12P23 + P12P33 + P22P13 + P22P33 + P32P13 + P32P23) 
+ 7 3 (P11P22P33 + P11P23P32 + P12P21P53 + P12P31P23 + P13P21P32 + P13P22P31)) ■ 



p = (Tr (Pi) + Tr (P2)) /°" 2 - The total signal to interference ratio is defined by ^ = Tr (Fx) /Tr (P2)- The 
spread of the signal power across the three BS locations is assumed to follow an exponential profile, as 
in G2ll . so that a range of possibilities can be covered with only one parameter. The exponential profile 
is defined by 

P tk = K k (a)a l ~\ (63) 

for receive location i £ {1, 2, 3} and source k where 

K k {a) = Tx{P k )/ (l + a + a 2 ), k = 1,2, (64) 

and a > is the parameter controlling the uniformity of the powers across the antennas. Note that as 
a the received power is dominant at the first location, as a becomes large (a 1) the third location 
is dominant and as a — > 1 there is an even spread, as in the standard microdiversity scenario. Using 
these parameters, eight scenarios are given in Table U for the case of two single antenna users. In Fig. [2l 
we explore the capacity of scenarios S1-S4 where ur = 3. The capacity result in (1371) agrees with the 
simulations for all scenarios, thus verifying the exact analysis. Furthermore, the approximation in (1571 ) 
is shown to be extremely accurate. All capacity results are extremely similar except for SI, where both 
sources have their dominant path at the first receive antenna. Here, the system is largely overloaded (two 
strong signals at a single antenna) and the performance is lower. The similarity of S3 and S4 is interesting 
as they represent very different systems. In S3, the two users are essentially separate with the dominant 
channels being at different antennas. In S4, both users have power equally spread over all antennas so the 
users are sharing all available channels. Figure [3] follows the same pattern with S6 (the overloaded case) 
being lower and the other scenarios almost equivalent. In Fig. |3j the overall capacity level is reduced in 
comparison to Fig. [2] as ^ = 10 implies a weaker second source. 

Figures |4] and [5] show results for a random drop scenario with M = ur = 3,W = N = 3 and n R = 



15 

6, M = 3, W = N = 6 respectively. In each random drop, uniform random locations are created for the 
users and lognormal shadow fading and path loss are considered where agF = 8dB (standard deviation 
of shadow fading) and 7 = 3.5 (path loss exponent). The transmit power of the sources is scaled so that 
all locations in the coverage area have a maximum received SNR greater than 3dB, at least 95% of the 
time. The maximum SNR is taken over the 3 BSs. Hence, each drop produces a different P matrix and 
independent channels are assumed. The excellent agreement between the approximation in (|57l) and the 
simulations in both Fig. |4] and Fig. [5]is encouraging as this demonstrates the accuracy of (1571 ) over different 
system sizes as well as over completely different sets of channel powers. Note that at high SNR, Fig. |5] 
gives much higher capacity values than Fig. |4] since there are 6 receive antenna rather than 3. In this high 
SNR region, the 7^ term in (|60l ) dominates and capacity can be approximated by log 2 (l + Perm (P) 7^) . 
With riR = N = 3 there are 6 cross products in Perm (P) whereas with tir = N = 6 there are 720 cross 
products. Hence, the bound clearly demonstrates the benefits of increased antenna numbers. In practice 
there is a trade-off between the costs of increased collaboration between possibly distant BSs and the 
resulting increase in system capacity. In Figs. [2] - [5l at low SNR the capacity is controlled by P T . Hence, 
since p = Pt/ct 2 , all four drops have similar performance at low SNR and diverge at higher SNR where 
the channel profiles affect the results. The upper bound and associated approximations are shown in Figs. 
|6] and |7J both for a two user scenario (S3) and a random drop. In Fig. |6j the upper bound is shown 
for scenario S3 as well as the high and low SNR approximations. The results clearly show the loss in 
accuracy resulting from the use of the simple Jensen bound. However, the bound is quite reasonable 
over the whole SNR range. The low SNR approximations are quite reasonable below dB and the high 
SNR version is as accurate as the bound above 15 dB. In Fig. |7J similar results are shown for a random 
drop with M = ur = 6, W = N = 6. Here, similar patterns are observed but the low and high SNR 
approximations become reasonable at more widely spread SNR values. For example, the low SNR results 
are accurate below dB and the high SNR results are poor until around 30 dB. In contrast, the upper 
bound is reasonable throughout. Hence, although there is a noticeable capacity error at high SNR, the 
cross-product coefficients in (|6T| ) and (|62l ) are seen to explain the large majority of the ergodic capacity 
behavior. 

VIII. Conclusion 

In this paper, we have studied the ergodic sum capacity of a Rayleigh fading macrodiversity MIMO- 
MAC. The results obtained are shown to be valid for both independent channels and correlated channels, 
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TABLE I 

Parameters for Figures[2]and[3] 





Decay Parameter 




Sc. No. 


User 1 


User 2 


? 


SI 


a = 0.1 


a = 0.1 


1 


S2 


a = 0.1 


a = 1 


1 


S3 


a = 0.1 


a = 10 


1 


S4 


a = 1 


a = 1 


1 


S5 


a = 0.1 


a = 0.1 


10 


S6 


a = 0.1 


a = 1 


10 


S7 


a = 0.1 


a = 10 


10 


S8 


a = 1 


a = 0.1 


10 



u 



Simulation 

Analytical approximation 

• Exact 




Fig. 2. Exact, approximated and simulated ergodic sum capacity in flat Rayleigh fading for scenarios S1-S4 with parameters: ur = 3, 
AT = W r = 2andc = l. 

which may occur when some of the distributed transmit/receive locations have closely spaced antennas. 
In particular, we derive exact results for the two-source scenario and approximate results for the general 
case. The approximations have a simple form and are shown to be very accurate over a wide range of 
channel powers. In addition, a simple upper bound is presented which demonstrates the importance of 
various channel power cross products in determining capacity. 



From (1331 ), lb can be written as 



Appendix A 
Derivation of I b 



dh 

001 



(65) 
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Fig. 3. Exact, approximated and simulated ergodic sum capacity in flat Rayleigh fading for scenarios S5-S8 with parameters: iir = 3, 
N = W = 2 and <; = 10. 
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Fig. 4. Approximated and simulated ergodic sum capacity in flat Rayleigh fading for M = tir — 3, W — N — 3 and four random drops. 



where 




-a 2 t-a 2 9 2 



>oJ U^Ai + tPiz + eiPuPa + esPii 



-d6 2 dt. 



(66) 



From (l66l) . L b becomes 




-a 2 t-a 2 9 2 



Jo U nR (O-, -I — ?2 r * 1 1_] 



■dOodt. 



(67) 
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Fig. 5. Approximated and simulated ergodic sum capacity in flat Rayleigh fading for tir = 6, M = 3, W = N — 6 and four random 
drops. 




Fig. 6. Ergodic sum capacity in flat Rayleigh fading for scenario S3 with parameters: M = tir — 3, W = N = 2 and ? = 1. 



Defining 



\P1P2\h 



(68) 



we use a partial fraction expansion in 9\ to give 



E 




A (6» 2 , e" 



■dQ 2 dt, 



(69) 
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Low SNR approx. 
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Fig. 7. Ergodic sum capacity in flat Rayleigh fading for a random drop with parameters: tir — 6, M — 3 and W — N — 6. 



where 



Ai (0 2 ,*) 

Pik 



1 



1 1 



Pk2 P'i2 
1 1 

-Pfei Pa 



Jik — Rk — Ri 

Ri 



Pi\Pi2 

To compute (l69l) . the following substitutions are employed 



(70a) 

(70b) 

(70c) 

(70d) 
(70e) 



u = a 2 t + a 2 9 2 
v . = _L + ^2_ 

P'il P%2 

The Jacobian of the transformation in (I71bl) can be calculated as 



Ji 



a 



1 1 

P%2 Pil 



Substituting (I71bl) and (O in © gives 

nR 




— ; -avi du , 

u Ji (vi + 9i + Ri) 



(71a) 
(71b) 



(72) 



(73) 
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where 



Ai (u, Vi) 



njji ( a ikVi + b ik U + J ik ) 



a 



a>ik — —r (ot ik — 0ik) 

b J_ f Pik_ _ Otik \ 
Ji \Pi2 Pil J 

The term Ai (u, Vi) in (I74al) can be written as a summation using partial fractions, to give 

B ik (u) 



A, t (u, = ^ 



' Vi + q ik u + r ik 



where 



B ik fit 



n R -3 



Cikl = hl^ik — CLilhk 

dikl = CLiklil — lik^il 

bik 

&ik 
Jik 

O-ik 



Qik 



Substituting (|75T) in (|73T) and simplifying gives 



n R n R r oo 




EE 

i=l ky^i 



First, we integrate over Vi in (1771) to give 



Bik (u) e~ u 



dvidu 



Ji (vi + 9i + Ri) (vi + q ik u + r ik ) 



n R n R poo 

'a 



EE 

i=l fc^i 



J; 



^ + 0i + i?i (A ifc M + r ik ) 



+9i + Ri\nikU + r ik ) 



where 



C ik {u, 6i) 

^ik 
Hik 



B ik 0) 



QikU + r ik - 0i - i?i 
1 

2 + Qik 



Pix<r 2 

1 



+ qik- 



(74a) 

(74b) 
(74c) 



(75) 



(76a) 

(76b) 
(76c) 
(76d) 
(76e) 



(77) 



(78) 



(79a) 
(79b) 
(79c) 



Let 



D ik (u, dx) = In 



[p^2 +0\ + Ri)(\ikU + r ik ) 
(1^2 +0i + R^\{iMkU + r ik ) 

then B ik [u] in (I76al) can be rewritten as the summation 



, k Cikiu + diki 
where 

/ \ n R— 3 

Substituting (EQ) and (l79cl) in <EHb gives 

«fl «fl njj „oo f , 

^ = EEE/ a*Mi)% — — u d \ 

^ ZT^ ,Tl -A; ( c ^« + «»fcj) + r ifc 

1=1 fe^l i^tljK 

Equation (l83l) can be further simplified to give 

" " Zm(M bikl -N h 



n R n R n R 

' ikl J ~" °ikl i 



L » = EEE 



/ J. A' / / /./.A- ' 

where 

D ik (u, 6{) du 



M 



b,. 



■ikl 



fl(°l) {U + Siki)' 

D ik (u,9x) du 



oc 



/l(fl) («+/2 («!))' 

and e iW = d ik i/c ik i. Next, we introduce the following linear functions of Q\. 

fl (#l) = ^ifei — C ik [9i 

Q 

h = m i« , 
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d 

Wi 



Dik MO 

h (0i) 



_ Cjkl 

j 1= o n ?fe; 



In 



ife« + 1 



1 



F«2 



dM h 



d0 1 



00 

Wi 



h (0i) 



du 



dN btkl 


r°° d 




du 




e 1= o Jo d °i 


[ fi(0i) . 


e 1= o ( U + m ikl) 



e 1= o ( u + e ^) 
Afc QM?i) 



ik 



e 1= o (u + mm) 



(91) 



(92) 



rdu (93) 



M h 



Cjkl 
n ikl 
b ikl 

niki 



Hi \ Ri , , 



P 2 a 2 



Hi(rik,eiki,Xik) - Hi [Ri,Siki, 

r 2 



Hi (rik,Siki,Hik) 



e pEEi\ — ) -e^Et (e iW ) 



£ ikl 

n-iki 



e^Ei[ — ) -e^Et (e iW ) 



(94) 



iV,, 



Cikl 
n iklHk 



Cjkl 
U ikl 



if 2 -Ri, Wjfc/, 



#2 (^fe , miki , Art ) - i? 2 -Ri , m lfei , 



H2 (fife, miki,fJ>ikj 



Hi Ri,miki, 



1 



iW 



+ Hi(rik,m lkl ,\ik) - Hi R4, m ik i , 



1 



Pi<x 2 



- Hi (r ik ,miki, fiik) 



+ 



m 



ikl 



niki 



e*ZE 1 [ — ) -e m ^Ei (m ik i) 



Uikl 



e n-2Ei[ — ) -e m ^Ei {m iH ) 



(95) 



where 



, Cikl 

niki — fikCiki — ctikiqik 5- 



miki 



(89) 
(90) 



bik QikRi 

Next, we can differentiate Mt ikl and Nb ikl and integrate over u to give the final result along with (l65l) and 
(1681) . Hence, from ([87]) and (I80al) we get (|9TT> . Substituting dED in d85]) and d86j) we get (HU) and (l93l) . 
(|92| ) and (|93T ) can be solved in closed form to give (|94| ) and (1931) , where we have used the two integrals 
defined as follows 



H\ (a, 6, c) 
H 2 (a, 6, c) 



00 e~* In (ct + a) 
Th6 
00 e~* In (ct + a) 



(t + by 



dt 



dt 



(96) 
(97) 
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and the constants are given by 



£ ikl 



m 



ikl 



£iki Pil 

1 



m ikl — p- 



b ikl 



> m ikl 



1 



1 



Both i?i and H 2 can be solved in closed form as 



H l (a, b, c) = e b E 1 (b)lnc + D 1 [~-b,b 



H 2 (a, b, c) = In c 



- - e b E, (6) 



2e 6 D x --6,6) + 



(f-6) L 



e% (6)-e?£ x - 
c 



where L>i (a, 6) is defined by 



£>i(a,6) 



00 e"* In (t + a) 



dt, for 6^0. 



Appendix B 



Calculation of E 



a 2 I 



Let Ai, A 2 , . . . , Afe_i be the ordered eigenvalues of H k H k . Since n R > (k — 1), all eigenvalues are non 
zero. Then, 



'fc-i 



E{\a 2 I + H k H k \}=E n( ff2 + A 

ffc— 1 



2\ fe— i— ll 



(98) 



where (|98T ) is from (flTT) and Lemma[2l Therefore, the building block of this expectation is E <j Tr^ ( i/ fe if fc 
From Lemma |2l 



Tr, (#ftf fc ) = £ 



^fc Hk 



(99) 



Therefore, from Lemma Q3 



i^Tr, (fffff fe )} = $>erm((Q fc r— 
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where the n R x (k — 1) matrix, Q;,, is given by 



E{k k oH k }=Q k . 
Note that summation in (1 1001) has (^y 1 ) terms. Then, the final expression becomes 

k— 1 

E{ I a 2 / + f if fc | } =J2Y: Perm ) (a 2 ) fc " i " 1 . 



(100) 



(101) 



i=0 a 



Appendix C 



a 



2 I + H k E A . 



Calculation of |S fc | £ 
A simple extension of (|45T) allows the expectation in the denominator of (1441) to be calculated as 



k-l 



where 



and from (fT4 



i=0 



^(t) = J]Perm((S fe U) ,Tl ' fe - 1 ) 



^ fe0 (*) = !■ 



The term in (11031) can be simplified using (fl"5T) to obtain 



Perm((Q fc )i fc - 1} 



(S fc ) 



Then, 



fc-i 



|EJ E 



i=0 



jfc-i-1 



where (£) = |E&| ^ (£). From (11041) . we obtain 



e«(t) = E|( E *)^J pem (w*)^ } 



(102) 



(103) 



(104) 



(105) 



(106) 
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where a nR -i nR is the compliment of a iyTlR . Therefore, it is apparent that £ ki (t) is a polynomial of degree 



n R — i. Clearly |S fc | E 



aH + H H k H k l H k 



is a polynomial of degree n R , since £ fc0 (t) = |S fc | is the 



highest degree polynomial term in t in (11051) . Then, 

i=o v 7 

Hence, applying (11071) in (11061) . 

riR—i / , \ I 



ii d — i ,n t 



(7 Z=0 v 7 



and ^ fc j (i) becomes 



riR-i / . \l 

z=o v 7 



where 



and from (fl4l) . ^ fc zo simplifies to give 

<^o = Tr z (P fc ) . 
Equation (11091) follows from (11081) due to the fact that 

Tr, ({P k ). . ) = for Z > n R - i. 
Therefore, (11021) can be written as 



which in turn can be given as 



ft— 1 n R 



2\ k—l—i—1 



i=0 1=0 



nR 

Z k \E{yi + H"x- k 1 H k \}=J2t l 

1=0 



(107) 



(108) 
(109) 
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where 

fc-i 

<pu = fan (v 2 ) 1 ■ (no) 

i=0 



Appendix D 
Extended Laplace Type Approximation 

Note the well-known fact that, a 2 I = E [A H A}, for an iid complex Gaussian matrix ensemble, A, of 
CAT ^0, random variables, where A is a k x k — 1 matrix as in [|22l . This result can be rewritten in 
the limit to give a 2 1 = K ^ a [A H A). Using this in (g3) gives 

A H A + H k H k 



h it) = — |— lim E 

\2-i k \ k— >oo 



A H A + H^H k 



1 



1 



lim E 



lim E 







(i) 













IS J K ~ y °° 



\ B k B k\ 

\Bj!T, k B k \ 



B ki] hi, 



where S fc = diag 2 j and B fc = V Using the well-known fact 

fc-i , 

from standard linear algebra, where b ki is the i* h column of B k , we can approximate (II 131) by 

/ / / II \ - I rr\ \ 



bki(l - B ki (B ki B k ^j B ki 



Bua \ bki 



l ~ l \b ki (l, k - T, k B ki (B ki E k B k ^ B ki T, k ^j b ki 



(HI) 



(112) 



(113) 



(114) 



(115) 



where b k i and B k correspond to a large but finite value of k. Approximation (II 151) assumes that the terms 
in the product in (II 14|) are independent. This is only true when b k i contains iid elements. However, in the 
macrodiversity case, all the elements of b k i are not iid. Nevertheless, part of b k i (the contribution from A) 
is iid. This motivates the approximation in (111 51 ) . Next, we apply the standard Laplace type approximation 
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11201 in (0131) to give 



fc -i J-B w [B u B ki \ B kl \b ki 



1 IV i 1 Kl 1 l*» V O nib J Kl I ~ M 

? g /_ — : : r i/ H _ x v (116) 

i=1 < b ki I Sfc — T, k B kl lB ki T, k B ki j B ki T, k J b k 







1 °ki 


(i - B ki (B ki Bk^j 


B k ^j fefcij 


> 


1^1 


nti 1 & ^ 




— ^>kB k i (B k JE, k B k i 


) B ki E k 





1 




\B k B k \ 


} 




E { \B k T, k B k 


1} 



(117) 



(118) 

Hence, a combination of approximate independence, the Laplace approximation for quadratic forms and 
the limiting version in (jl 1 1|) gives rise to the approximation used in Sec. |V] The accuracy of this approach 
is numerically established in the simulation results in Sec. IVIII 
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